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Abstract 

Insoluble surfactant monolayers at the air/water interface undergo a phase tran- 
sition from a high-temperature homogeneous state to a low-temperature demixed 
state, where dilute and dense phases coexist. Alternatively, the transition from a di- 
lute phase to a dense one may be induced by compressing the monolayer at constant 
temperature. We consider the case where the insoluble surfactant monolayer inter- 
acts with a semi-dilute polymer solution solubilized in the water subphase. The 
phase diagrams of the mixed surfactant/polymer system are investigated within 
the framework of mean field theory. The polymer enhances the fluctuations of the 
monolayer and induces an upward shift of the critical temperature. The critical 
concentration is increased if the monomers are more attracted (or at least less re- 
pelled) by the surfactant molecules than by the bare water/air interface. In the 
case where the monomers are repelled by the bare interface but attracted by the 
surfactant molecules (or vice versa), the phase diagram may have a triple point. 
The location of the polymer special transition line appears to have a big effect on 
the phase diagram of the surfactant monolayer. 
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1 Introduction 



Understanding the subtle interactions between macromolecules, such as polymer or pro- 
teins, and amphiphiles, such as surfactants or phospholipids, has been a problem of prime 
interest in recent years in many industrial applications and in biological systems. For 
instance, biomembranes |], § are usually depicted as fluid bilayers composed of different 
constituents: phospholipids, cholesterol and proteins. In addition, a complex macro- 
molecular network (the cytoskeleton) is associated with the inner side of the bilayer and 
modifies the mechanical properties of the membrane, while the glycocalix, on the outer 
side, is believed to play an important role in molecular recognition. In industry, surfac- 
tants are used in a wide range of applications (detergents, soaps, oil recovery, paints), 
where polymers are often added in order to provide stability for the system, especially in 
the case of colloidal suspensions and oil/water emulsions @. Those mixed polymer and 
surfactant systems tend to create complex self-assembly structures (connected micelles, 
gels, networks, etc.) f|]-f|. Finally, drug delivery via micro-encapsulation is an other ex- 
ample where the stability of surfactant vesicles is improved by the adsorption of polymer 
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In recent years, a new category of associating polymers has been introduced. Those are 
the hydrophobically modified water soluble polymers (HM-WSP), consisting of a water 
soluble polymer backbone carrying small hydrophobic side chains. Such polymers present 
interesting properties of self-association, which may even be enhanced by the addition 
of surfactant, and are very useful as viscosity modifiers of aqueous solutions || §]. The 
subtle coupling between the surfactant and the polymer may lead to unusual phenomena 
like thermogelation ||, where gelation of the system is obtained upon increasing of the 
temperature. Such systems have been studied theoretically as well as experimentally 
in the bulk. However, little is known about their behavior at interfaces [[II] . 



In the present work, the interaction of water soluble polymers with a surfactant mono- 
layer located at the air/water interface is considered. We restrict ourselves to the relatively 
simple situation of an insoluble surfactant monomolecular layer (Langmuir monolayer). 
Langmuir monolayers have been used in many applications [0, [13], from evaporation 
control to non-linear optic devices (via the creation of Langmuir-Blodgett monolayers). 
They are also used to study crystallization of solids [14j] and provide useful model systems 
for more complicated fluctuating liquid interfaces (membranes) where curvature effects 
can not be neglected. 

Another motivation for the present study comes from the lack of understanding of 
adsorption (or depletion) of polymers close to non-ideal interfaces, as compared with 
adsorption on ideal surfaces (namely, perfectly flat and chemically homogeneous). On 
ideal surfaces [p!5||-||19|| theories for neutral and flexible polymers in good solvent have 



been performed (both for adsorption and depletion) and compared with scaling theories. 
For non-ideal surfaces, much less theoretical works exists. It has been suggested that 
the bending properties of a curved interface are modified when in presence of adsorbing 
polymer [p0|1 -|p2|. When the polymer adsorbs on both sides of the interface (a bilayer 



for instance), the curvature modulus decreases, while the saddle-splay modulus increases. 
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When it adsorbs only on one side, a non-zero spontaneous curvature is induced. The 
situation of a perfectly flat but chemically heterogeneous interface has been considered 
only in a few works ||22||— 1|25||. The case of annealed disorder (namely, when the disorder 
is at thermodynamic equilibrium and the heterogeneities can diffuse laterally) is found to 
behave differently from the case of quenched disorder (where heterogeneities are spatially 
"frozen"). However, in both cases the adsorption of the polymer is increased by the non- 
ideality of the surface. In this context, a surfactant monolayer is an example of a non-ideal 
annealed surface, where the order parameter is the local surface surfactant concentration. 

The phase diagram of surfactant monolayers can be constructed as a function of the 
thermodynamical variables |l2| : surface pressure and temperature (or equivalently area 
per molecule and temperature). At low surface pressure, a phase separation occurs (for 
temperatures below the corresponding critical temperature). Dilute (gaseous) and dense 
(liquid-expanded) regions of the monolayer coexist, in analogy to phase transitions in 
the bulk. In the phase diagram, single-phase and two-phase regions are separated by a 
coexistence curve. At higher surface pressure, other phase transitions occur. Depending 
on the symmetries of the specific surfactant molecules, the phase diagrams are more 
complex and still a topic of current investigation |27J . 

In the following, we consider how a simple condensation transition (gas to liquid 
expanded) of a surfactant monolayer at the air/water interface is affected by the presence 
of polymer in the water subphase. The free energy and the assumptions used in deriving 
it are introduced in Sec. 2, while in Sec. 3 we discuss the main results, as applied to a 
simple case; the general theory is detailed in Appendices A and B. Finally, some analytical 
considerations on the critical point are presented in Appendix C. 



2 The polymer/surfactant free energy 

The model used for the mixed surfactant /polymer system follows closely the lattice model 



introduced in Ref . [24 . The local (dimensionless) free energy per site F, rescaled in units 
of fcsT, is the Boltzmann constant and T is the temperature, can be separated into 
three parts: the surfactant contribution F s , the polymer contribution F p , and the coupling 
term F ps : 

F = F S + F p + F ps (2.1) 
In the following, those three terms are discussed separately. 



2.1 The surfactant contribution F s 

The monolayer free energy is calculated using a lattice-gas model. Each lattice site is 
occupied either by a surfactant molecule or by an artificial vacancy, in order to allow us 
to consider a compressible monolayer. The free energy of a surfactant monolayer is the 
sum of the enthalpy and entropy of mixing and depends on the monolayer area fraction 
(or equivalently coverage) c ranging from zero to one, c = A /A, where A is the close- 
packing area of a surfactant molecule (or the area of one site on the lattice) and A is 
the actual area per surfactant molecule on the interface. Typically Aq ~ 25 — 35 A 2 for 
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a surfactant molecule Disregarding linear terms, the surfactant free energy F s (per 
site and per ksT), within a Bragg- Williams (mean field) theory, is written as: 

F s = v- l c{l - c) + clog c+ (1 -c)log(l-c) (2.2) 

where v~ x is the dimensionless interaction parameter of the surfactant on the surface 
and describes Van der Waals interactions between neighboring particles. The interactions 
between the head groups of the surfactant molecules, playing an important role in the 
determination of the highly compressed phases, as well as the freedom of conformation for 
the hydrophobic chains, whose coupling with the surfactant average c is determinant in 
the liquid expanded versus liquid condensed transition p8fl , are not taken into account. As 
only short-ranged interactions are considered (between neighboring sites), the surfactant 
molecules are supposed to be neutral. The main interaction modeled by the parameter v 
are the Van der Waals interactions. 

For an insoluble monolayer, the total number of surfactant molecules is fixed. At 
low (and positive) values of v (corresponding to low temperatures), a phase separation 
between dense and dilute regions follows from eq. ( |2.2j ). The stability of such a monolayer 
is obtained by studying the convexity of the free energy |[29|| . The monolayer becomes 
unstable if the second derivative of the free energy becomes negative. The condition 
F"(c) = defines the spinodal line, separating metastable and unstable regions. The 
spinodal line obtained from eq. (|2.2|) is v®(c) = 2c(l — c) and it lies within the coexistence 
region of the phase diagram. In addition, the coexisting curve limiting the two phase 
region (the binodal line) is easily found from eq. ( |2.2|) as the system is symmetric about 
c = 0.5: 

logc — log(l — c) 

The spinodal and binodal lines join together at the critical point c = 0.5, v c = 0.5. In 
Fig. 1 the binodal line and the critical point are shown for a pure surfactant monolayer. 

2.2 The polymer contribution F p 

The polymer in the sub-phase is assumed to be neutral and flexible as well as in good 
solvent conditions, hence with a positive second virial coefficient and no polymer-solvent 
phase separation. For a semi-dilute polymer solution, a mean field theory applied to the 



Edwards density functional method is commonly used |3(J, |19|| . The free energy density 
is conveniently expressed as a function of the variable 4>(z) related to c p (z), the local 
monomer concentration, by <f) 2 (z) = c p (z)/q,. The coordinate z denotes the perpendicular 
distance from the interface, and q, = c p (z — > oo) is the concentration of the polymer in 
the bulk (acting as a reservoir). The characteristic length in the solution is the Edwards 
correlation length £ = a/ ^3vCb, where v is the excluded volume parameter (positive, in 
good solvent conditions), and the typical energy parameter (per ksT) for the interactions 
between monomers in the bulk is e p = A ^v c\. Using these notations, the free energy per 
site reads: 

F P = ^fdzU(V<P) 2 + -^ 2 -l) 2 ) (2-4) 
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The first term accounts for the elastic flexibility of the polymer chains and the second 
originates from the excluded volume interaction combined with the equilibrium condition 
with the polymer bulk reservoir. The polymer free energy F p is a functional of the polymer 
profile <f)(z) and of the order parameter at the interface cj) s = <fi(z = 0). It does not include 
the energy of interaction with the surface, discussed separately below. 

Minimizing F p with respect to the polymer profile 4>(z), leaving the surface value as 
a free parameter, yields the polymer profile <j>{z) = coth(z/£ + b), where b is a constant 
of integration related to <p s by <p s = cothfr, in the case of adsorption (<f)(z) > 1), and 
<p(z) = tanh(2:/£ + b') in the case of depletion (<f)(z) < 1), where, similarly, <p s = tanhfr'. 
For both adsorption and depletion, the free energy F p for the optimal profile is 

F P= e f {fi - 30 s + 2) = | (<f> s - l) 2 (</>. + 2) (2.5) 

and has a minimum at 4> s — 1. This means that the polymer solution would like to be 
homogeneous throughout the solution at the imposed bulk value c p (z) = c&. The only 
possibility of obtaining a profile with <p s ^ 1 is due to the short range coupling of the 
polymer with the surface. This coupling includes the surfactant monolayer as well as the 
bare air/water interface. It is given below by the term F ps . 

A quantity accessible to experiment which measures the total adsorption of the monomers 
at the interface is the polymer surface excess defined as Y = dz(c p (z) — c&). Using the 
above results of the minimization for the polymer profile (mean- field theory), it is simply 
related to <p s by 

T = c^{(j> 3 - 1) (2.6) 

Note that, if p is the volume fraction of the monomers in the bulk solution, a naive 
calculation starting from c& ~ p/a? (where a is the size of a monomer) yields t p ~ p 3 / 2 . 
For a semi-dilute polymer solution, 7V -4//5 <C p 1 (where N is the number of monomers 
in a chain). Hence, roughly, for iV = 10 4 the range for typical values of e p is given |yj by 
1(T 4 < e p < 10" 1 . 

Although the self-consistent field theory provides a convenient and qualitatively correct 
framework for the description of the semi-dilute polymer solution, some of its predictions 
(like the form of the polymer profile c p (z) for instance) are in disagreement with a scaling 
theory |I9[ . Nevertheless, we will use it to model the polymer behavior in solution. 

2.3 The coupling term F ps 

A bilinear term in the surfactant and monomer concentrations at the interface (z = 0) is a 
simple, yet meaningful, phenomenological coupling for the polymer-interface interaction, 
which is assumed to be short-ranged: 

Fp, = -\[aoc + 7 o(l - cM = -^e ps (c - c*)0 2 (2.7) 

«o is the polymer/surfactant interaction parameter, and 70 is the polymer/bare interface 
interaction parameter. In eq. (|2.7|) , we define the " effective" polymer/surfactant inter- 
action parameter e ps = 0:0 — 70- It is positive whenever the monomers interact more 
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favorably with the surfactant molecules than with the bare water/air interface. The spe- 
cial transition coverage is defined as c* = — 7o/(«o ~ 7o)- In principle, those parameters 
can depend on temperature. 

The phenomenological coupling F ps can be justified for polymers in the semi-dilute 
regime since in such systems the monomers concentration is small with respect to unity. 
However, it represents only the lowest term in an expansion in the surfactant concentration 
at the interface. When e ps > 0, the F ps term corresponds to a repulsion of the polymer 
from the surface (depletion) for c < c* and to an attraction to the surface (adsorption) 
for c > c*. The special transition line c = c* occurs for physical values of the coverage, 
< c* < 1, when the attraction (repulsion) of the monomers with the surfactant molecules 
is in competition with the repulsion (attraction) with the bare interface. In the < c* < 1 
range, a positive «o is equivalent to having a positive e ps and means that the interaction 
between the monomers and the surfactant molecules is attractive. 

2.4 The total free energy F 

Combining all three contributions, F s + F p + F ps , we obtain the total free energy (per site 
of the interface and per ksT): 

F = z/ _1 c(l - c) + clogc+ (1 - c)log(l - c) 

+ l -t p {<f>, - i f (0 S + 2) - l -t ps {c - c*)(f)l (2.8) 

Note that the energy is invariant under the transformation e ps — *> —e ps , c — > 1 — c and 
c* — > 1 — c*. Therefore, it will be assumed in the following that e ps > without loss of 
generality. 

The free energy F is a function of <p s and c. Minimizing it first with respect to the 
polymer surface order parameter <p s (mean field approximation), we obtain 



2T p + \ 




(2.9) 



where A = e ps (c — c*) measures the strength of the interaction between the polymer and 
the overall interface (including the bare interface as well as the surfactant). Equation 
( OP relates 4%(c), the concentration of monomers at the interface, with the surfactant 
area fraction c. Consequently, the entire polymer profile and the polymer surface excess 
can be found as a function of the surfactant concentration on the interface. 

The limit of a very strong adsorption, A/e p — >• oo (e.g., c* 0), yields (f) s ~ X/e p 1. 
On the other hand, the limit of very strong depletion, X/e p — > — oo (e.g., c* ^> 1), yields 
(p s ~ \e p /\\ <^ 1. For \/e p — > (e.g., c — > c*), 4> s — > 1, which means that the polymer 
solution remains homogeneous: c p (z) = c&, even at the surface. As is shown in Fig. 2a, 
the special transition line c = c* divides the parameter range into an adsorption region 
(c > c*) and a depletion one (c < c*). 

If the surfactant monolayer is in the two-phase region, dilute and condensed regions 
of the surfactant coexist, and the polymer adsorbs differently on those regions because 
of its different affinity as described by the parameter e ps > 0. Note that as the curve 
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4> s (c) is convex (see Fig. 2b), the polymer surface excess is enhanced when the surfactant 
monolayer undergoes a phase separation. Qualitatively, the convexity of the curve s (c) 
indicates |24| that the surfactant concentration fluctuations increase the average polymer 
surface concentration S and, hence, the polymer surface excess r ~ (fi s — 1. 

Equation ( [2.9|) shows how the surfactant molecules affect the polymer adsorption. The 
main remaining task is to understand how the polymer itself affects the phase diagram of 
the surfactant monolayer, since the two problems are coupled. Using eq. Q2.9|) , (fi s = </> s (c), 



the total free energy can be written only as a function of the surfactant coverage c: 

F(c) = F s (c)-|(<^ + 3)0 s (2.10) 
o 

The study of the convexity of the total free energy F(c) as a function of c (the only re- 
maining order parameter) determines the location of the modified spinodal line. Similarly, 
the full phase diagram can be obtained numerically from a common-tangent construction 
of F{c). 

For sake of clarity, the main physical features can be illustrated in a simple situation. 
This is done in the next section, where e p , c* and ve ps are taken to be independent of the 
temperature T, and v ~ T. A more general treatment without any assumptions on the 
T-dependence of the interaction parameters is presented in Appendices A and B. 



3 Results 

3.1 T-dependence of the interaction parameters 

T jv is independent of the temperature if we assume that the interaction potential of the 
surfactant molecules has an infinite repulsive core followed by a weak attraction inde- 
pendent of the temperature, < w <C ksT. In such a case, the surfactant second virial 
coefficient is given by v s = 1 — w/ {ksT). Expanding the energy F s , eq. (|2.2j), in powers of 
c shows that the virial expansion |29| is equivalent to the Bragg- Williams theory (at low 
concentrations) with z/ _1 = w/(A;^T). Thus ksT/u is independent of the temperature. 
Note that the strength of this attraction w is related to the critical temperature of the 
pure surfactant monolayer by w = v~ x ksT c = 2/cbT c . 

For the semi-dilute polymer solution, taking only excluded volume interactions be- 
tween monomers and assuming an athermal solvent, v and e p are independent of T. 

We also assume that the interaction between the monomers and the bare interface 
is independent of the temperature, resulting in 7 ~ 1/T since the dimensionless 70 
is rescaled in units of ksT. In a similar manner, neglecting the steric effects between 
the monomers and the surfactant molecules, and assuming that the monomer /surfactant 
interactions are attractive, weak and short-ranged, results in a ~ 1/T. Under these 
conditions ve ps and c* are independent of T. 

In conclusion, the phase diagram in the simplified case can be plotted in the (z/, c) 
plane and is a cut of the global phase diagram (presented in the appendices) plotted 
in the (z/, c, e p , ve vs , c*) space. The next subsections will present some features of this 
simplified case. 
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3.2 The (is, c) phase diagram 

We limit ourselves to v > 0, e p > and e ps > 0; v > corresponds to an attraction 
between the surfactant molecules; e p > follows from the assumption that the polymer is 
in a good solvent (its excluded volume parameter v is positive), and e ps > can be used 
without loss of generality, as was explained in Sec. 2.4. The spinodal line v s (c) of the 
mixed surfactant-polymer system is obtained from the condition 

d# ~ + C(1-C) £ p #(c) + l " U V'V 

The critical point is the extremum point on the spinodal, satisfying in addition 

d 3 F , , 

a? = ° < 3 - 2 > 

Equation ( |3 . 1| ) shows that the spinodal temperature v s is shifted upwards [fM] and the 
region of instability is increased. Physically, this general effect comes from the indirect 
attractive interaction between the surfactant molecules induced by the polymer, and was 
already explained elsewhere fl2"2 |. Here, it is represented by the term — e p (0 3 + 30 s )/6 in the 
free energy. It is bigger for larger values of the surfactant concentration c because </> s (c) 
is an increasing function of c. Consequently, the phase diagram is no longer symmetric 
around c = 1/2. The shift on the spinodal and binodal lines is bigger for the large values 
of c and the critical point is shifted to a concentration c c > 1/2. 

In the following, we discuss several limits and try to show that (except in the low- 
coupling limit) the position of the special transition line c = c* in the (u, c) plane crucially 
affects the phase diagram. This can be checked easily in the limits of very high and 
very low surfactant concentration, where the scaling of the spinodal temperature v s is 
analytically derived and depends on the relative position of the spinodal line with respect 
to the special transition line. 

3.3 The limits c -> and c 1 

Three cases can be distinguished in the c — > limit of the spinodal line: 

• if c* > 0, the polymer is strongly depleted from the interface. The shift of the 
spinodal temperature v s {c) from the pure value z^!?(c) is small and scales as v s — ~ 
c 3 , where = 2c(l — c) is the pure surfactant spinodal line (no added polymer). 

• On the other hand, if c* < 0, the polymer is strongly attracted by the interface and 
v s scales as v s ~ c 1//3 . 

• In the special case when c* = 0, the polymer solution remains homogeneous. The 
dominant term in the spinodal equation also changes and v s scales as u s ~ c 1//2 . 

The scaling of the spinodal temperature in the limit c — > 1 depends similarly on the 
position of the special transition line relatively to the line c = 1. In the following, the 
other regions of the phase diagram are considered and v is at least of order unity. 
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3.4 The low-coupling limit: ue ps <C 1 



In the case when ve ps is small enough ]52 , the third term in the spinodal equation ( |3.1| ) 
is negligible. The phase diagram is very similar to that of the pure surfactant monolayer 
(see Fig. 1). The spinodal temperature v s (c) as well as the binodal temperature z/ b (c) can 
be expanded in powers of ve ps . To second order in ve ps , the shift of the binodal line is 
identical to the shift of the spinodal line: 

v s {c) - i£(c) = v b {c) - u° b (c) = ^{^ P sY (3.3) 

For the low-coupling limit, the critical temperature is shifted upwards [[24[] by a factor of 
order (^e ps ) 2 . Note that, to second order in ve ps , the shift is independent on the surfactant 
area fraction and the special transition coverage c*. The critical concentration c c is also 
shifted upwards, but only to third order in z/e ps , : 

c c -l/2 = -^(ve ps ) 3 (3.4) 



Another case where the third term in the spinodal equation ( |3.1| ) is negligible is when 
4> s (c) <C 1. Here, the polymer is very strongly depleted from the interface. This occurs, 



for instance, when e p C 1 and c < c* 



3.5 The strong coupling limit: large values of ue 



ps 



We define the strong coupling limit as the situation when the spinodal line is strongly 
shifted upwards, v s (c) ^> 1 and the first term in the spinodal equation ( |3.1| ) is negligible. 
e ps is the only parameter left in the equation depending on the temperature. As v = 
(ye ps )/e psi the temperature on the spinodal line is proportional to the parameter ue ps 



and the coefficient of proportionality e 1 is obtained from the spinodal equation ( |3.1| ) 
rewritten as a fifth order equation for e ps depending on e p , c* and c (and independent of 





C C* 3 4 (c C*) 3 - 



-f A - f 4 - v > r 5 = n 

^■ps ^-ps J) \i -i „\ ps \ J 



c{l-c) ps ps e 2 p c(l-c) ps 

Whenever eq. ( |3.5| ) has a unique positive solution e ps , the approximation of the strong 
coupling limit is self-consistent provided that ue ps ^> e ps . This is the case for c > c*. 
On the other hand, when the polymer is depleted from the interface (corresponding to 
c < c*), there is a minimal value of c for which eq. ( |3.5|) has a positive solution. Indeed, 
a sufficient condition for eq. ( |3.5| ) not to have any positive solution is: 

( 11 V /4 
c* - c > = JT P ~ 0.2524,/i; (3.6) 

V 1350 + 210v / 42/ 



This implies that, in a system where the condition fl3.6p is respected for the largest physical 
c value, c = 1, the strong coupling limit can not be defined and the critical temperature 
v c is necessarily of order unity. 
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When c* — > — oo or when e p — > (and c > c*), the solution of eq. ( J3.5p is small, 
e ps < 1. In this situation, i/e ps of order unity is enough to ensure that v s 1 and: 

v.{c) = ue ps (c - c*) l ^[c(l - c)] 1 /^; 2 / 3 (3.7) 

In the strong coupling limit, the binodal line is also proportional to the coupling 
parameter ve ps . Numerical solution of eqs. (|3.1| )-( j3T2"D indicates that the behavior of the 
binodal line is simpler than the one of the spinodal line. In particular, if the strong 
coupling limit can be defined only over a range of surfactant coverage c (for instance, 
when the special transition line intersects the phase diagram), then the binodal line is 
found to be proportional to ve ps for all values of c (except for the limiting values c — > 
or for c — > 1). Figure 3 shows the dependence of the binodal temperature on the special 
transition value c*, at a fixed value of the polymer interaction parameter e p . Lower values 
of the special transition coverage c* correspond, for a fixed coupling ve pS) to a higher 
binodal temperature. 

3.6 High polymer flexibility: e p <C 1; the possibility of a triple 
point. 

As was explained in Sec. 2, e p is likely to be very small for a semi-dilute polymer solution. 
When e p <C 1, the strong coupling limit is obtained for c > c* (as was explained in Sec 
3.5) and the polymer is strongly adsorbed at the interface (see Fig. 3). On the other 
hand, for c < c*, the polymer is depleted from the interface and the situation is a one of 
low coupling. Consequently, the features of the phase diagram are particularly sensitive 
to the position of the special transition line: when c* > 1, the phase diagram is very close 
to the one of the pure surfactant monolayer (see Fig. 1), while for c* < the spinodal line 
is given by eq. ( |3.7| ) and the binodal temperature is strongly shifted upwards. When the 
special transition line intersects the phase diagram, there is a competition between the 
two types of phase behavior. In Sec. 3.5, it was mentioned that the phase diagram of the 
strong coupling is dominant in that situation. 

More precisely, an expansion of the free energy in powers of e p /e ps shows that, for c > 
c*, F(c) ~ — l/Q(e 3 ps /e 2 p )(c — c*) 3 , while for c < c*, it yields F(c) ~ F s (c). Fig. 4a shows a 
plot of the free energy in a case where this sudden decrease of F(c) for c > c* is particularly 
well defined. When c* > 1/2, the binodal line of the pure surfactant monolayer can be 
obtained with a common tangent construction (in the range of temperatures where the two 
coexistence concentrations are lower than the special transition coverage c*). However, it 
corresponds to a metastable state. There is yet a second pair of coexisting concentrations 
which leads even to a lower free energy. 

When < c* < 1, the phase diagram of the strong coupling limit predominates over 
the one of the pure surfactant monolayer, while for c* > 1 this is not the case anymore. 
For physical values of the parameters, the transition between the two regimes, occurring 
around c* = 1, is smooth. In some situations, the two pairs of coexisting concentrations 
mentioned above are stable and correspond to two first-order phase transitions, as is shown 
in Fig. 4b. The total phase diagram exhibits two critical points and one triple point. An 
example of such a phase diagram with a triple point is presented in Fig. 5. In Fig. 6 
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we show for the same set of parameters three typical isotherms plotted in the reduced 
surface pressure II - reduced area per surfactant plane. The reduced surface pressure 
is defined by rescaling the actual surface pressure (i.e., the difference between the bare 
water/air surface tension and that of the surfactant monolayer) by t4 ^//cbT, resulting 
in the relation II = uc 2 d(F/c)/dc. Depending on the temperature, the isotherms can 
have zero, one or two plateaus corresponding respectively to zero, one or two coexisting 
regions. 



4 Discussion 

In order to derive our model we assumed several simplifying assumptions. The expression 
for the surfactant free energy F s does not take into account the surfactant hydrophobic 
tails degrees of freedom. The coupling between the tail conformations and the concentra- 
tion of the surfactant molecules becomes crucial at high surfactant densities and can lead 



to a rich phase behavior |26|, [27]], which was not addressed here. We consider here only 
the dilute phases of surfactant monolayers: the gaseous phase at low densities and tem- 
peratures, the liquid-expanded phase at higher densities and the condensation transition 
between them ||28|| . 

Another limitation of the model comes from our mean-field treatment of the polymer 
free energy F p . The assumption that the polymer solution is semi-dilute may break 
down close to the interface if the polymer strongly adsorbs. Our approach assumes that 
c p (z = 0) <C 1 or equivalently A < p, where p is the volume fraction of the monomers 
in the bulk and A = e ps (c — c*). The coupling term F ps should be regarded as the first 
term in an expansion. A more precise study should take into account higher order terms, 
particularly in the surfactant concentration c. 

For water soluble polymers, hydrogen bonds play an important role because of the 
strong polarity of water molecules and more refined expressions for F p and F ps should be 
used. For example, the phase diagram of water soluble polymers like PEO exhibits closed 
loops of immiscibility and the definition of good solvent conditions is somewhat vague || . 

We are not aware of many experiments performed on adsorption of polymers on a 
Langmuir monolayer, which will allow a direct comparison with our results 0. Inter- 
esting polymers that can be used experimentally to test our theory are hydrophobically 
modified water soluble polymers (HM-WSP) [ffl,|j8fl-|ll||. These associating copolymers 
are attracted by the bare interface because of short hydrophobic side chains, attached 
covalently to the main chain. By adjusting the number and the length of the side chains, 
one can directly modify their surface affinity In our model it corresponds to 70, the 
parameter of interaction between the polymer and the bare surface. 

Figure 7 shows two different kinds of interaction between the monomers and the in- 
terface. In (a), monomers interact with the interface through any kind of short-ranged 
interaction: either attractive or repulsive. In (b), the attraction of HM-WSP polymers 
towards the interface is illustrated: one of the aliphatic side chains of a HM-WSP polymer 
is in a low energy state in the air subphase. The monomers of the main chain covalently 
bound to this aliphatic group are consequently attracted in the region of the interface. 

The HM-WSP polymers can be used to study systematically the dependence of the 
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phase diagram on the special transition concentration c* (as well as on the coupling param- 
eter e ps ), including the interesting situation where the special transition line intersects the 
phase diagram (0 < c* < 1) because the polymer is repelled from the surfactant molecules 
(«o < 0). The surfactant can be chosen as nonionic, with a polar head (hydrophilic) iden- 
tical to the monomers, resulting in a repulsion between the "brush" (formed by the polar 
heads) and the polymer (water being a good solvent for both) [F2^| . The resulting coupling 
parameter e ps is negative and the triple point occurs for values of c* close to zero and not 
close to unity as when e ps > (Fig. 5). However, such polymers usually have a compli- 
cated behavior in the bulk and their self-assembly properties may be crucial. Another 
possible candidate for experiments may be a polymer with a "simpler" behavior in water, 
like PEO, having a surface affinity |33]]. Probably here one needs to treat more explicitly 
the hydrogen bonds. 

The parameters in the model can also be tuned by changing the surfactant character- 
istics. Pentadecanoic acid, for example, whose gaseous to liquid expanded transition has 
already been studied by several authors can probably be an interesting surfactant 

to use. 

In a previous work [fm , several phase diagrams for the mixed polymer-surfactant 
system have been proposed from qualitative and general arguments. The proposed phase 
diagrams exhibit the special transition line for the polymer (c = c*) and the coexistence 
line for the surfactant molecules, but the interaction between these two lines was not 
treated in detail. Our results suggest that the position of the special transition line has 
a very strong effect on the position of the coexistence line. Ref. also predicted a G 
line separating a region where the polymer segregates from the surfactant from a region 
where the polymer and surfactant are mixed together. Our mean-field model does not 
address directly this prediction since we have only one minimum of the polymer surface 
value 4> s as a function of surfactant coverage c [see eq. (|2.9|)], and we assume that the 
polymer solution is homogeneous in the directions parallel to the interface as long as the 
surfactant monolayer is homogeneous. 

iFrom our study it seems that the O line of Ref. p2] and the coexistence line are 
the same: regions of different concentrations for the polymer correspond to regions of 
different concentrations for the surfactant. The polymer can be attracted to the interface 
and segregate from the surfactant only if it is attracted by the bare interface but repelled 
by the surfactant molecules. This situation is driven by energy terms which are first 
order in the polymer concentration [see eq. (|2.7|) 1 and is consequently different from the 



one predicted in Ref. P5 |. It will be interesting to understand in a detailed way this 
discrepancy, especially for < c* < 1, where the coupling between the special transition 
and the surfactant phase diagram leads to the richest variety of phenomena. 



5 Conclusion 

We addressed in detail the adsorption of a semi-dilute polymer solution on a surfactant 
monolayer, and the resulting phase diagrams. In our model, the most important degree 
of freedom is the local concentration c of surfactant at the interface. Since the monomers 
interact with regions of different concentrations with an energy proportional to c — c* (c* 



12 



is the concentration at the special transition) , a rich variety of phenomena results from the 
coupling between the polymer solution and the surfactant monolayer. The polymer surface 
excess is enhanced and the phase diagrams of the surfactant monolayer are modified. 
The monomers induce an additional indirect attraction between the surfactant molecules 
depending on the concentration of surfactant on the interface. Consequently, the range 
of the homogeneous region in the phase diagram decreases. When the monomers interact 
more favorably with the surfactant molecules than with the bare interface, the critical 
concentration itself increases. 

The value of the special transition coverage c*, describing the interaction of the 
monomers with the surfactant molecules relatively to their interaction with the bare wa- 
ter/air interface, has a major effect on the phase diagram of the surfactant monolayer. 
When the monomers are repelled by the surfactant molecules as well as by the bare in- 
terface (for instance, e ps > and c* > 1), the phase diagram is not very different from 
the one of a simple surfactant monolayer without polymer in the water subphase. On the 
other hand, when the monomers are attracted by the surfactant molecules as well as by 
the bare interface (for instance, e ps > and c* < 0), the increase of the two-phase region 
can be important. In the intermediary situation both when the polymer is attracted to 
the interface or depleted (0 < c* < 1), these two scenarios are in competition leading to 
a rich phase behavior and in some cases, the phase diagram displays two critical points 
and one triple point (as in Fig. 5 and 9k). 

Finally, we mention two possible extensions of the present study. The first is to consider 
the situation of soluble surfactants, and to take into account the complex surfactant- 
polymer interactions in the bulk. Interesting experimental results have been obtained 



for such sy stems PHI . Moreover, our model can easily be adapted to mixed monolayers 
with two components. Such binary mixtures have been shown to exhibit immiscibility 
at room temperature, (e.g., for cholesterol and DMPC PBL B7I]). The understanding of 



the features of this phase transition may help to understand interactions of polymers 
with flexible membranes (lipid bilayer) and the phenomenon of budding of membranes 
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A Appendix A: The Spinodal Equation 



The general equation of the free energy of the system, as was derived in Sec. 2, is 

F = v- x c{\ - c) + clogc + (1 - c) log (1 - c) - ^ (fi + 3<j) s - 4) (A.l) 

If we do not assume any specific temperature dependence for z/, e p , e ps and c*, the global 
phase diagram has to be studied in the five dimensional space of z/, e p , e ps , c* and c. 
The binodal as well as the spinodal surfaces are four- dimensional hypersurfaces in this 
five-dimensional space. We shall consider both positive and negative v. It is worthwhile 
to investigate the spinodal surface because it can be done analytically, and it roughly 
describes the phase separation region, since the spinodal surface lies always within the 
coexistence region. The surface of instability (spinodal surface) also indicates the locations 
of possible critical points. For fixed values of e p , e ps and c*, the spinodal equation ( |3.1| ) 
describes cuts of the spinodal surface by the curve v" 1 = ^ _1 (c) 

It is interesting to look also at other cuts of the parameter space beside the (z/ _1 ,c) 
direction. An analytical expression can also be found in the (e p , c) direction (but not for 
the other interaction parameters). 
Defining u as 

u = °-ZSL (— L_ - 2uA (A.3) 
e ps \c{l — c) J 

we note that u = for both the special transition line and the spinodal line of the non- 
coupling case (F(c) = F 8 (c)). From the spinodal equation ( |A.2|) it follows that, if e p is a 
solution of the spinodal equation ( |3.1|) , e 2 is the solution of a bi-quadratic equation: 

Au 2 e 4 p - A (l + iu - u 2 ) e\ - X A u = (A.4) 



Alternatively, it can be shown that, if e p is a solution of ( A.4 ), then, either e p or — t p 
is a solution of the spinodal equation (ft.lQ. As the polymer is supposed to be in good 
solvent conditions, we disregard negative solutions of the spinodal equation for e p . The 



second order equation for e p ( |A.4 ) may have either only one real solution, or even no real 



solution at all, depending on its discriminant. Furthermore, only positive solutions should 
be taken into account, since ei > 0. All these considerations give several bounds to the 
possible values of u and A. The roots of the bi-quadratic equation (|A.4j ) are: 



\2 

f3 ± = —,(l + Au-u 2 ±(u + l)Vl + 6u + u 2 ) (A.5) 



When P + is real and positive, only one solution of the spinodal equation exists: either 
e p = y/m. or e p = — y/J3+. We shall consider this solution only in the case when it is 
positive. The same can be said for 

Defining u = — 3 + v8 ~ —0.1716, we separate three cases for the roots of ( A.4 ) given 
by eq. flO) : 
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• For u > 0, f3 + is real and positive (while /3_ is negative, and therefore, irrelevant). 
Hence, there is only one solution for e p . 

• For uo < u < 0, f3 + and /3_ are both real and positive (for u — > 0, /3_ — > 0). Hence, 
there are two solutions for e p . 

• For -u < Mo, neither (3+ nor /3_ are real and positive. Hence, no physical solution for 

dp. 

Once e ps and c* are fixed, e p is a function of v and c. It is instructive to identify the 
domains in the (z/, c) plane where the spinodal equation (|3.1|) has zero, one or two real 
solutions. This requires to locate the lines u = and u = u . Depending on the values 
of e ps and c*, two different situations are identified and shown in Fig. 8. In particular, it 
can be seen in Fig. 8b that, for 1/2 < c* < 1, there can be a range of values for v~ l where 
cuts of the spinodal surface in the plane (e p , c) are disconnected (if e ps is below a certain 
critical value). For instance, in Fig. 9j (1/2 < c* = 0.86 < 1), a cut of the phase diagram 
at Uq 1 < u^ 1 = 1.8 < is presented. The spinodal line is defined for two disconnected 
regions of concentration: one centered around c = 0.6 and one at high concentrations. 

B Appendix B: The Global Phase Diagram 

The phase diagram for the surfactant monolayer, including the surfaces of instability and 
coexistence, is presented in the three dimensional parameter space of z/ _1 , e p and c through 
several cuts in two dimensions (e p and c), while fixing the other two parameters, e ps and 
c*. We first discuss some general features of the phase diagram, starting with the simple 
case of no effective coupling between the polymer and the interface (e ps = 0; F(c) = F s (c)) 
as a reference point. The dependence of the phase diagram on the position of the special 
transition line is then studied, and the phase diagrams are presented for the cases: c* < 0, 
c* > 1 and < c* < 1. 

When there is no effective coupling between the polymer solution and the surfactant 
monolayer, e ps = and the special transition line is not defined. The concentration of the 
polymer at the interface is identical to its bulk value and depends on the strength of the 
interaction between the monomers and the interface = 70 = A. The phase separation 
in the monolayer was explained in Sec. 2.1, and all its features are independent on the 
polymer parameter e p . In the plane (z/ _1 ,c), the line v~ x = l/[2c(l — c)] delimits the 
unstable region (see Fig. 8). From eq. ( |A.2| ) it appears that as soon as there is a non zero 
effective coupling, the line of instability occurs for smaller values of v . It means that 
the coupling enlarges the region of instability of the homogeneous state. Consequently, 
the dark shaded region inside the line z/ _1 = l/[2c(l — c)] in Fig. 8 is always a zone of 
instability, for any value of the parameters e p , e ps and c*. When e p is big, the polymer 
"stiffness" induces the polymers in solution to be uniform, and the coupling acts as a 
chemical potential having a value \e pS) as can be seen from eq. ( |2.10| ) (after expanding cf) s 
to first order in powers of e^ 1 ). In this limit, the properties of the phase separation are 
the same as in the non-coupling case. Considering the cuts of the phase diagram in the 
direction (e p , c) for e ps = 0, there are two different regions of the parameter space: 
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• the region v~ x < 2, where there is no instability and no phase separation. 

• the region u^ 1 > 2, where there is a phase separation. 

Fig. 9a shows the spinodal and the binodal lines for v = 2.31. Both lines are vertical 
because the properties of the phase transition are independent of e p , as was already 
explained above. In the limit v^ 1 — > 2, the phase diagram reduces to a line of critical 
points parallel to the e p axis at c = 0.5. The coupling results in a deformation of the 
line of critical points and enlarges the phase coexistence region. The cuts of the phase 
diagram in the direction (e p , c) for z/ -1 < 2 are topologically very different from the cuts 
obtained for v^ 1 > 2, because the latter necessarily contain the region of instability of 
the non-coupling situation (even for large values of e p ) , while, for the former, large values 
of e p necessarily correspond to domains of stability of the monolayer. 

We first discuss the cases where the special transition line c = c* does not intersect 
the parameter space (c* < or c* > 1). The phase diagrams are constructed by numerical 
solution of eqs. Q-Q. 

B.l The & < case 

When both the surfactant and the interface attract the monomers, two situations can 
occur: if ao > 7o > 0, then e ps > and c* < 0. But if 70 > «o > 0, e ps < and 
c* > 1. However, as was explained at the beginning of Sec. 2.4, these two situations 
can be mapped onto each other. Therefore, we concentrate here on the case e ps > and 
c* < 0. The only relevant solution of the spinodal equation for e p is then e p = +y//3+ , 
defined outside of the region of instability of the non-coupling case. When e p is big, the 
polymer stiffness enforces the polymer solution to be uniform, S ~ 1, and the monomers 
exert a strong uniform attraction on the surfactant molecules. 

For z/ _1 < 2, there is a phase separation, provided that e p is small enough. This is 
shown in Fig. 9b, where it can be seen that the spinodal and the binodal lines join at 
a critical point. When u~ l is changed, the basic topology of the binodal and spinodal 
lines remains unchanged. But, in the limit v^ 1 — > 2 _ , close to the region where a phase 
separation occurs even without a coupling, the value of e p at the critical point increases. 
Moreover, the critical concentration tends to the value of the non-coupling case c = 1/2 
(because e p at the critical point is large). If, on the other hand v~ x decreases, the value of 
e p at the critical point decreases and the critical concentration c c increases, c c — > 1 when 
v v — ► —00. 

For z/ -1 > 2, the phase diagram is presented in Fig. 9c. There is no critical point 
because the phase separation occurs for all values of e p . For e p — > 0, the effects of the cou- 
pling between the monomers and the surfactant molecules are the strongest: the polymer 
is strongly attracted to the interface, particularly in the regions dense in surfactant. As 
a result, the surfactant molecules aggregate with a close-packing coverage c = 1. 

B.2 The c* > 1 case 

The opposite situation happens when both the surfactant and the bare interface repel 
the monomers. As explained above it is enough to consider only e ps > and c* > 1. In 
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this situation, the function ^ _1 (c) defined by the line u = uq in the plane (i/ -1 ,c) (see 
Appendix A) has a minimum v~ x < 2, which is a function of e ps and c*. When z/ -1 < z/^ 1 , 
then -u < Mo- The spinodal equation has no solution and the homogeneous monolayer 
state is stable. For z/ -1 > i/" 1 , from a topological point of view, the spinodal surface 
is deformed for intermediate values of e p . Whereas for large values of e p , the polymer 
remains stiff and the properties of the phase transition are those of the non-coupling case, 
for small values of e p , the repulsion of the interface is dominant. The polymer is strongly 
depleted and the coupling has no effect. 

The cuts of the phase diagram in the direction (e p , c) have the following features: 

• For Vq 1 < u^ 1 < 2 the phase diagram has a closed loop immiscibility curve with an 
upper and a lower critical points as shown on Fig. 9d. When v~ x — > 2~, the upper 
critical point tends to oo and the lower one to 0, whereas when v~ x — > i/" 1 , the 
domain of instability becomes very small. 

• For z/ -1 > 2, there is a phase separation, whose characteristics differ from the one 
of the non-coupling case only for intermediate values of e p . An example of such a 
case is presented in Fig. 9e. 

The line of critical points has only, for intermediate values of e p , a small distortion 
with respect to the straight line of the non-coupling case. One consequence is that there 
is a maximal value for the critical concentration. When c* is increased, the repulsion from 
the interface increases and the phase diagram resembles the non-coupling case, even at 
intermediate values of e p . 

B.3 The < & < 1 case 

We turn now to the more complex case where the special transition line intersects the 
physical range of parameter space (0 < c < 1). The situation for < c* < 1/2 is rather 
simple, and two typical phase diagrams are shown in Figs. 9f (for e ps < 2) and 9g and 
(for e ps > 2). They resemble the c* < case. However, for 1/2 < c* < 1, depending 
on the value of the coupling parameter e ps , the competition between the different terms 
in the free energy leads to a large variety of possible phase diagrams. An example of 
a phase diagram with a reentrant phase is presented in Fig. 9h. In Fig. 9i yet another 
type of a phase diagram is presented where a second metastable region resides inside the 
two phase region. In some cases, as was suggested at the end of Appendix A, the phase 
diagram appears as if it is composed of two disconnected parts (those are cuts obtained 
for Uq 1 < z/ -1 < 2). An example is given in Fig. 9j. Note that for 2 < z/ -1 < z/f 1 the 
spinodal surface is disconnected while the binodal surface is connected, since at very low 
values of e p , a phase separation always occurs between very dense (c ~ 1) and very dilute 
(c ~ 0) regions (for c* < 1), as is explained below. In the three dimensions (e s , e p , c), the 
phase diagram is always connected: when z/ _1 is increased the disconnected parts of the 
phase diagram meet and this may result in the appearance of a triple point for a range of 
values for e s (an example is shown in Fig. 9k). These triple points occur only for a small 
range of values for e p ; typically, for e p ~ 0.1. 

Note that in all these phase diagrams (Figs. 9f to 9k), the binodal line at low values 
of e p always ends at c = and c — 1. It can be understood since for low values of e p , 
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regions of surfactant concentrations c < c* (and outside of the region of instability of the 
non-coupling case) are not unstable because they repel the monomers, and the coupling 
consequently does not have any effect. The free energy is unchanged as compared to the 
non-coupling case. However, it is strongly modified for c > c*, as has been explained 
in Sec. 3.6. and is shown on Fig. 4a. The common tangent construction results in a 
couple of stable coexistence points (c ~ 0, c ~ 1). Therefore, the regions of surfactant 
concentrations c < c* (and outside of the region of instability of the non-coupling case) 
are metastable. 



C Appendix C: The Critical Point 

It is possible to obtain analytical results for the critical point, by solving simultaneously 
eqs. ( |3.1|) and ( p. 2D , for the two unknowns: the critical concentration c c and the critical 
temperature. 

C.l Bounds on the critical concentration 

In general, no analytical expression for the critical concentration c c exists but we derive 
an analytical expression for an upper bound. Equation (|3.2j ) can be rewritten as a third 
order equation for (j) 2 s : 

# + 3# + 3/(c)# + /(c) = (C.l) 
where the coefficient /(c) is defined as: 

1 — 2c 

/(C) = {l-2c) + c\l-cfelJel (C ' 2) 

We denote C\ as the value of the surfactant concentration for which f(c) diverges. A 
careful study of eq. flU.lp shows that /(c) < 0, and yields a non- analytical solution for (fi s 
and consequently for c c . It can also be shown that 1/2 < c c < c\. From its definition, c\ 
is the solution of a fourth order equation and has a (complicated) analytical expression 
as an odd increasing function of eL/eS. For e 3 /e 2 — 0, ci = 1/2 and when eL/e? — > oo, 



ci->i m 



C.2 The critical point in the strong coupling limit 

Hereafter we concentrate on the specific case of the strong coupling limit, where v~ x can 
be neglected in the spinodal equation eq. 



c(l 

e 



b\ + 1 = e 2 p J 3 s ^ (C.3) 



Taking the polymer order parameter e p and the special transition c* as known parameters, 
the characteristics of the critical point can be determined by solving a system of three 
equations (the spinodal equation (|C.3|) , eq. (|3.2j ) and the definition of 4> s ( |2.9|) ) with three 



unknowns: the concentration of surfactant c, the temperature appearing through e ps and 



20 



<f) s (which has been added for mathematical convenience). First, bounds on the critical 
concentration are obtained, depending only on the special transition concentration. Then, 
a method of determination of the critical concentration is discussed. Using eq. (|C.3|) , 
eq. ( |3~2| ) can be rewritten as a second order equation for 2 . 

As there is at least one real and positive solution to this equation, a bound on the 
value of the critical concentration, depending only on the special transition concentration 
c* is obtained (Fig. 10). 



For & < 0: 1 - 
For < c* < 1: 



Vc* 2 - c* + 1 <c c <\ (c* + 1) + 



*2 



1 



(c* + 1) - Vc* 2 - c* + 1 



<c c <| 



(c* + 1) + Vtf 



.★2 



c* + 1 



For c* > 1: 



;c* + 1) - Vc* 2 - c* 



i 



< c c < 1 - c* + Vc 



*2 



C* + 1 



For c* = 0, this bound is < c c < 2/3. Therefore, if the effective coupling parameter e ps 
is positive, we know that the critical concentration obeys 1/2 < c c < 2/3. 

The second order equation for 2 can be rewritten as a second order equation for the 
critical concentration c c , depending on the special transition concentration c* and on the 
concentration of monomers at the interface d> q : 



[4 



3(0 2 + l) 2 ]c 2 -2 



c*{<fc + I) 2 







(C.4) 



The study of this equation shows that: 

• For c* = or c* = 1, it reduces to a first order equation. 

• For c* < 0, c c = c_ = 
(2-(l + c*)(0 2 + l) 2 - 

• For c* > 1, c c = c + = 



l) 4 (c* 2 



c^ + l)-4(0 2 + l) 2 + 4 /(4-3(0 2 + l) 2 ) 



(2 - (1 + c*)(0 2 + l) 2 + ^(02 + 1)4(c , 2 _ c* + i) _ 4(02 + 1)2 + 4 ) / (4 _ 3(0 2 + 1)2) 

For < c* < 1, the situation is more complicated. It can be shown, in particular, 
that there is a minimum for the concentration of the monomers at the interface at 
the critical value. This minimal value always corresponds to a situation of depletion 
for the polymer solution. It can also be shown that in some situations, when e ps > 
and 1/2 < c* < 1 or symmetrically when e ps < and < c* < 1/2, the solutions c„ 
and c + may be relevant: this accounts for the possibility of two critical points. 



The definition of 



H) can be rewritten as: 

,9 ^ 



-ps 



[C — C 



(C.5) 



e ps can be eliminated from the spinodal equation ( |C.3|) by using its expression obtained 
from eq. (U.5). On the other hand, the concentration can be substituted by using its 
expression c_(0 s ,c*) or c + (0 s , c*). One is left with a high order polynomial equation for 
4> s depending only on the special transition concentration c* and the polymer interaction 
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parameter e p . Once this equation is numerically solved, the other characteristics of the 
critical point (concentration, temperature) are easily obtained as simple functions of <j) s . 
In particular, once the critical concentration is known, the critical e ps , hence the critical 
temperature are obtained from eq. (C.5). 



22 



Figure Captions 



Figure 1: The phase diagram for a bare surfactant monolayer, without polymer in 
the water subphase. At low v (low temperatures), the homogeneous state is unstable 
and the binodal line delimits the two-phase coexistence region labeled A+B. The 
critical point is located at v c = 0.5, c c = 0.5 and is shown by a dot. 

Figure 2: (a) The polymer order parameter <p s and (b) the polymer surface excess, 
T = r/(c&£) = 4> s — 1, as function of the surfactant coverage c for e ps /e p = 10 and 
c* = 0.5. In (a) the special transition line divides the region where the polymer 
is adsorbed (c > c*) from the region where it is depleted (c < c*). In (b), for 
homogeneous monolayers of concentrations C\ = 0.29, C2 = 0.875 and c* = 0.5, 
the polymer surface excess is respectively —0.6, 3 and 0. Due to the convexity of 
the curve, for a surfactant monolayer of average concentration c* demixing between 
coexisting regions of concentration c\ and C2, the surface excess T is positive. 

Figure 3: The value of e~} on the binodal line in the strong coupling limit (inde- 
pendent on v) for c* = (dashed line), c* = 0.5 (dotted line) and c* = 0.8 (full line) 
at a fixed value of the polymer interaction parameter e p = 0.1. The binodal temper- 
ature v is directly related to e p } by v = {ye ps )e~}. When c* decreases, the polymer 
interacts more favorably with the interface and its effect on the phase diagram of 
the surfactant monolayer is stronger: the binodal temperature increases. 

Figure 4 : The free energy F as a function of the surfactant concentration c at the 
interface for v = 0.48, vt ps = 1, c* = 0.95. In (a) e p = 0.02 and in (b) e p = 0.1. The 
common tangent construction is shown by thin lines and the coexistence values by 
circles. In (a), in the region < c < c*, the plot is almost similar to the case of the 
pure surfactant monolayer. The corresponding coexistence region, 0.33 < c < 0.67 is 
only metastable, since it is contained in the second coexistence region 0.11 < c < 1.0. 
In (b), both coexistence regions 0.31 < c < 0.72 and 0.83 < c < 0.99 occur. 

Figure 5: The phase diagram of the surfactant monolayer for c* = 0.95, e p = 0.1 
and ve ps = 1. The two-phase region labeled A+B ends at the critical point: v c = 
0.51, c c = 0.50. The second B+C critical point is located at v c = 0.53, c c = 0.94. 
All three two-phase regions: A+B, A+C and B+C join at a triple point: v = 0.46, 
C4 = 0.25, c B = 0.79 and c c = 0.99, where all three phases (A, B and C) coexist. 
Critical points are shown by a dot. 

Figure 6: Isotherms for the surfactant monolayer. The reduced surface pressure 
II = uc 2 d(F/c)/dc is plotted versus the reduced area per molecule l/c—1 = A/A — l 
on a logarithmic scale. Aq is the close-packing area and F is the total free energy 
as defined in the text. Three typical isotherms are shown for three different tem- 
peratures: v = 0.56 (dotted line, no phase transition), v = 0.52 (dashed line, one 
phase transition), ^ = 0.49 (full line, two phase transitions). The other parameters 
are identical to the ones of Fig. 5. 



Figure 7: Two different kinds of interaction between the monomers and the inter- 
face. In (a), monomers interact with surfactant molecules either attractively (for 
instance, as a consequence of favorable van der Waals interactions) or repulsively (if, 
on the contrary, water is a good solvent for both molecules, and the surfactant heads 
act as a polymer brush). In (b), an associating polymer is shown; the hydrophobic 
sub-chains tend to dispose themselves in the air subphase and consequently attract 
the polymer chain close to the interface, (adapted from Ref. [^2| ) 

Figure 8: Regions of definition of the spinodal hypersurface in the plane (— z/ -1 , c) 
for: (a) e ps = 1 and c* = 0.8; (b) e ps = 1 and c* = 0.95. The full line is the 
u = uq line defined in Appendix A. The curved dashed line is the line of instability 
for the non-coupling case: 2z/ _1 = l/[c(l — c)]; the vertical dashed line is the line 
c = c*. The region where only the solution f} + is real and positive, and generates a 
positive solution e p to the spinodal equation is light shaded, while the region where 
both solutions (3 + and /3_ are real and positive, and generate a positive solution e p , 
is hashed. The dark shaded region is the zone where the surfactant monolayer is 
always unstable. The region above the full line (u = uq) is a zone where neither j3 + 
nor /?_ are real and positive. In (b), the u = u$ line has two extrema at v — v and 

v — V\. 

Figure 9: Cuts of the phase diagram in the (e p , c) plane. The full lines represent 
the binodal surface, while the dashed lines represent the spinodal surface. Critical 
points are shown by a full circle, (a) z/ _1 = 2.31 and e ps = 0; this plot shows 
the reference situation of the non-coupling case. Next plots illustrate the text of 
Appendix B, in the three cases c* < 0, c* > 1 and < c* < 1. In (b) to (j) 
e ps = 1. (b) is- 1 = 1 and c* = -1; (c) v~ l = 2.2 and c* = -1; (d) u~ l = 1.95 and 
c* = 2. This plot shows a closed loop phase diagram, (e) v~ x = 2.2 and c* = 2. 
The e p scale in this plot is logarithmic for convenience, (f) v~ x = 1 and c* = 0.3; 
(g) v^ 1 = 2.2 and c* = 0.3; (h) u^ 1 = 1.8 and c* = 0.84. This phase diagram 
displays a reentrant phase, (i) u^ 1 = 2.1 and c* = 0.84. The e p scale in this plot 
is logarithmic for convenience and the shaded region corresponds to the metastable 
states, (j) u~ l = 1.8 and c* = 0.86. This plot shows a disconnected phase diagram. 
In (k) v~ l = 2.1, e ps = 2.1 and c* = 0.95. A first two-phase region A+B extends to 
infinite values for e p ; the second two-phase region ends at a critical point c c = 0.90, 
i e p)c = 0.16; all three two-phase regions: A+B, A+C and B+C join at a triple point: 
e p = 0.06, ca = 0.30, cb = 0.70 and cc = 0.99, where all three phases (A, B and C) 
coexist. 

Figure 10: In the strong coupling limit, the range of the critical concentration 
c c depends only on the special transition concentration c*, and is shown inside the 
shaded region. 
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